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intervals for the parameter of the Poisson distribution is discussed. 
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1 Introduction 

In paper the unified approach to the construction of confidence intervals 
and confidence limits for a signal in the background presence, in particular, 
for Poisson distributions, is proposed. The method is widely used for the 
presentation of physical results Q though a number of investigators criticize 
this approach |^]. 

Here we use a simple method of constructing confidence intervals for the 
Poisson distribution parameter for a signal in the presence of background 
which has the Poisson distribution with the known value of parameter to 
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compare with a conventional procedure Q . The method is based on the state- 
ment that the probabihty of true value of the Poisson distribution parameter 
to be a specified value (in the case of the observed number of events x) dis- 
tributes in accordance with a Gamma distribution. It is shown that this 
statement has both Bayesian and frequentist interpretations. The experi- 
mental results often give non-integer values for a number of observed events 
X (for example, after the background subtraction Q) when the Poisson dis- 
tribution occurs. That is why there is a necessity to have a procedure for 
constructing the confidence intervals in this case. The paper offers a gener- 
alization of Poisson distribution for a continuous case ^ . The generalization 
given here allows one to construct confidence intervals and confidence lim- 
its for the Poisson distribution parameter (for integer and real values of a 
number of observed events) using conventional methods. 

In Sect. |2| the interrelation between the frequentist and Bayesian def- 
initions of the confidence interval is shown. The method of constructing 
confidence intervals for the Poisson distribution parameter for a signal in 
the presence of background which has the Poisson distribution with the 
known value of parameter is described in Sect. ^ The results of confidence 
intervals construction and their comparison with the results of the unified 
approach are also given in Sect. |^. In Sect. ^ the generalization of Poisson 
distribution for the continuous case is introduced. The examples of confi- 
dence intervals construction for the parameter of the Poisson distribution 
analogue and for the Poisson distribution parameter using the Gamma dis- 
tribution are considered in Sect. ^ and in Sect. |6|. The main results of the 
paper are summarized in the Conclusion. 

2 The interrelation between frequentist and Bayesian 
definitions of confidence interval 

Let us have a random value taking values from the set of numbers x € X. 
Consider the two-dimensional function 

f{x,X) = ^e-\ (2.1) 



^The early version of the study can be found in S.I. Bityukov, N.V. Krasnikov, 
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where x > and A > 0. 

Assume, that the set X includes only the whole numbers, then for each 
value of A a discrete function /(x, A) describes the distribution of probabili- 
ties for the Poisson distribution with the parameter A and a random variable 
x, i.e. ^ ~ Pois{X). 

Let us write down the density of Gamma distribution Ta^^+i as 

/(x,a,A) = -^^e-''^A^ (2.2) 
T[x + 1) 

where a is a scale parameter, x + 1 > is a shape parameter, A > is a 
random variable, and T{x + 1) is a Gamma function. Since the x is integer, 
then x\ = r(x + 1). Note that this notation is also used in the case of real 
X. Let us set a = 1, then for each x a continuous function 

/(x. A) = ^e"^, A > 0, x> -I (2.3) 

is the density of Gamma distribution Ti^^+i with the scale parameter a = 1 
(see Fig. ||). The mean, mode, and variance of this distribution are given by 
x + 1, x, and X + 1, respectively. 

Assume that in the experiment with a fixed integral luminosity (i.e. the 
process under study is considered as a homogeneous process for a given time) 
the X events of a Poisson process are observed. It means that we have an 
experimental estimation A(x) of the parameter A of the Poisson distribution. 
We have to construct a confidence interval (Ai(x), A2(x)), covering the true 
value of the parameter A of the distribution under study with a confidence 
level 1 — a, where a is a significance level. It is known from the theory 
of statistics [^, that the mean value of a sample of data is an unbiased 
estimation of the mean of distribution under study. In our case the sample 
consists of one |^ observation x. 

For the discrete Poisson distribution the mean coincides with the esti- 
mation of parameter value, i.e. A = x in our case. 

Let us consider the formula 

P{X\x) = P{x\X) = ^e~\ (2.4) 
x! 

^The Poisson distributed random values have a property: if ^ ^ Pois{Xi) and r; ~ 
Pois{\2) then ^ + ~ Pois{Xi + A2). It means that if we have two measurements £1 and 
X2 of the same random value ^ ~ Pozs(A), we can consider these measurements as one 
measurement xi + X2 of the random value 2 ■ ^ ~ Pois{2 ■ A). 
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This formula (2.4) results from the Bayesian formula |^ 

P{X\x)P{x) = P{x\X)P{X) (2.5) 

in the assumption that all the possible values of parameter A have equal 
probability, i.e. -P(A) = const. In this assumption the probability that 
unknown parameter A obeys the inequalities Ai < A < A2 is given by the 
evident Bayesian formula 

P(Ai < A < A2IX) = P(Ai < X\x) - P{\2 < X\x) = P{X\x)dX, (2.6) 

where P{X\x) is determined by formula (2.4). 

Formula (2.6) has also a well defined frequentist meaning. Using the 
identity 

°° \*p-Ai /.A2 \£p-A ^ \^p-^2 

i=x+l -^^1 i=0 

one can rewrite formula (2.6) as 



^'(Ai < A < X2\x) = \-Pin < x\X2)-P{n > x\Xi) = P{n < x\Xi)-P{n < x\X2), 

(2.8) 

where P{n < x\X) = ^ j — and P{n > x\X) = ^ j — . 

n=x+l 

The right hand side of formula (2.8) has a well defined frequentist mean- 
ing and it is the definition of the confidence interval in the frequentist ap- 
proach. Note, that this definition of the confidence interval for the Poisson 
distribution parameter is self-consistent both for the case Ai = A2 and for 
the case x = 0. As an example of the shortest 90% CL confidence interval of 
such type in case of the observed number of events x = 4 is shown in Fig. 

In ref. Q (pp. 406-407) the interrelation between the frequentist and 
Bayesian definitions of confidence interval is shown, nevertheless, the author 
criticizes the Bayesian approach of the confidence interval determination. 

As it is seen from the identity (2.7) the probability of true value of 
parameter of Poisson distribution to be equal to the value of A in the case of 
one measurement x has probability density of Gamma distribution Fi^i+f. 

^Notice that this identity takes place for any Ai > and A2 > (in particular, if 
A2 < Ai). 
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Correspondingly, in the case of m measurements xi,X2, ■ ■ ■ , Xm of the 
random values ■^i, ^2, ^m; where ~ Pois{X) for i = 1,2, ... ,m, the 
probability of true value of parameter of Poisson distribution to be equal to 
the value of A has probability density of Gamma distribution i+J^"" xi ■ 

3 The method of confidence intervals construction 

Let us consider the Poisson distribution with two components: Signal com- 
ponent with a parameter and background component with a parameter 
Af,, where A;, is known. To construct confidence intervals for the parameter 
As of a signal in the case of observed value x, we must find the distribution 

P{\s\x). 

Firstly let us consider the simplest case x = s + h = 1. Here s is the 
number of signal events and h is the number of background events among 
the observed x events. 

The h can be equal to and 1. We know that the h is equal to with 
probability 

PO = P{b = 0) = ^^e-^'' = e-^" (3.1) 
and the h is equal to 1 with probability 

pi = P{h = 1) = ^e-^" = X^e-^K (3.2) 

Correspondingly, P(6 = Olx = 1) = Pis = l\x = 1) = — — — and 

Po+Pi 

P{h = l\x = 1) = P{s = Olx = 1) = — . 

It means that the distribution of P{Xs\x = 1) is equal to the sum of 
distributions 



p{s = i\x = i)ri,2+P(s = o\x = i)ri,i = ^^ri,2+^^ri,i, (3.3) 

Po+Pi Po + Pi 

where Fi^i is the Gamma distribution with the probability density ^(Asls = 
0) = e"^" and Fi^2 is the Gamma distribution with the probability density 
P(A5|s = 1) = A^e"^^. As a result, we have 

P(A,|x = l) = ^^e-^^ (3.4) 
J- + Ah 
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Table 1: 90% C.L. intervals for the Poisson signal mean As, for total events 
observed x, for known mean background Af, ranging from to 2. Comparison 
between the results of ref.[l] and the results from the present paper. 



x\Xb 


0.0 ref.[l] 


0.0 


1.0 ref.[l] 


1.0 


2.0 ref.[l] 


2.0 





0.00, 2.44 


0.00, 2.30 


0.00, 1.61 


0.00, 2.30 


0.00, 1.26 


0.00, 2.30 


1 


0.11, 4.36 


0.09, 3.93 


0.00, 3.36 


0.00, 3.27 


0.00, 2.53 


0.00, 3.00 


2 


0.53, 5.91 


0.44, 5.48 


0.00, 4.91 


0.00, 4.44 


0.00, 3.91 


0.00, 3.88 


3 


1.10, 7.42 


0.93, 6.94 


0.10, 6.42 


0.00, 5.71 


0.00, 5.42 


0.00, 4.93 


4 


1.47, 8.60 


1.51, 8.36 


0.74, 7.60 


0.51, 7.29 


0.00, 6.60 


0.00, 6.09 


5 


1.84, 9.99 


2.12, 9.71 


1.25, 8.99 


1.15, 8.73 


0.43, 7.99 


0.20, 7.47 


6 


2.21,11.47 


2.78,11.05 


1.61,10.47 


1.79,10.07 


1.08, 9.47 


0.83, 9.01 


7 


3.56,12.53 


3.47,12.38 


2.56,11.53 


2.47,11.38 


1.59,10.53 


1.49,10.37 


8 


3.96,13.99 


4.16,13.65 


2.96,12.99 


3.18,12.68 


2.14,11.99 


2.20,11.69 


9 


4.36,15.30 


4.91,14.95 


3.36,14.30 


3.91,13.96 


2.53,13.30 


2.90,12.94 


10 


5.50,16.50 


5.64,16.21 


4.50,15.50 


4.66,15.22 


3.50,14.50 


3.66,14.22 


20 


13.55,28.52 


13.50,28.33 


12.55,27.52 


12.53,27.34 


11.55,26.52 


11.53,26.34 



Using formula (3.4) for P{Xs\x = 1) and formula (2.8), we construct the 
shortest confidence interval of any confidence level in a trivial way. 

In this manner we can construct the distribution of P{Xs\x) for any 
values of x and A;,. As a result, we have obtained the known formula |9| 

P(A.|x) = ^^4^e-^^ (3.5) 
i=o 

The numerical results for the confidence intervals and the results of pa- 
per [Q] are compared in Table H and Table ||. 

It should be noted that in our approach the dependence of the confidence 
intervals width for the parameter A^ on the value of Xb in the case x = is 
absent. For x = the method proposed in ref. [^0| also gives a 90% upper 
limit independent of A;,. This dependence is also absent in the Bayesian 
approach 11|. 
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Table 2: 90% C.L. intervals for the Poisson signal mean A^, for total events 
observed x, for known mean background A^ ranging from 6 to 15. Compar- 
ison between the results of ref.[l] and the results from the present paper. 



x\Xb 


6.0 ref.[l] 




6.0 




12.0 ref.[l] 




12.0 




15.0 ref.[l] 


15.0 





0.00, 0.97 





00, 2 


30 


0.00, 0.92 





00, 2 


30 


0.00, 0.92 


0.00, 2.30 


1 


0.00, 1.14 





00, 2 


63 


0.00, 1.00 





00, 2 


48 


0.00, 0.98 


0.00, 2.45 


2 


0.00, 1.57 





00, 3 


01 


0.00, 1.09 





00, 2 


68 


0.00, 1.05 


0.00, 2.61 


3 


0.00, 2.14 





00, 3 


48 


0.00, 1.21 





00, 2 


91 


0.00, 1.14 


0.00, 2.78 


4 


0.00, 2.83 





00, 4 


04 


0.00, 1.37 





00, 3 


16 


0.00, 1.24 


0.00, 2.98 


5 


0.00, 4.07 





00, 4 


71 


0.00, 1.58 





00, 3 


46 


0.00, 1.32 


0.00, 3.20 


6 


0.00, 5.47 





00, 5 


49 


0.00, 1.86 





00, 3 


80 


0.00, 1.47 


0.00, 3.46 


7 


0.00, 6.53 





00, 6 


38 


0.00, 2.23 





00, 4 


19 


0.00, 1.69 


0.00, 3.74 


8 


0.00, 7.99 





00, 7 


35 


0.00, 2.83 





00, 4 


64 


0.00, 1.95 


0.00, 4.06 


9 


0.00, 9.30 





00, 8 


41 


0.00, 3.93 





00, 5 


15 


0.00, 2.45 


0.00, 4.42 


10 


0.22,10.50 





02, 9 


53 


0.00, 4.71 





00, 5 


73 


0.00, 3.00 


0.00, 4.83 


20 


7.55,22.52 


7 


53,22 


34 


2.23,16.52 


1 


70,16 


08 


0.00,13.52 


0.00,12.31 



4 The Generalization of Discrete Poisson Distri- 
bution for the Continuous Case 

Let us consider the case when x € X are the real values and denote x\ = 
r(x + 1), then we can consider the function 

/(x,A) = ^e-^ (4.1) 

as a continuous two-dimensional function. Fig. ^ shows the surface described 
by this function. Smooth behaviour of this function along x and A (see Fig. ^) 
allows one to assume that there is such a function /(A) > — 1, that 

/ f{x,X)dx = l (4.2) 

Jl{X) 

for the given value of A. It means that in this way we introduce a contin- 
ued analogue of Poisson distribution with the probability density f{x, A) = 
A^ 

— j-e~^ over the area of the function definition, i.e. for x > 1{X) and A > 0. 
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The values of the function /(x, A) for integer x coincide with corre- 
sponding magnitudes in the probabiUties distribution of discrete Poisson 
distribution. Dependences of the values of function Z(A), the means and the 
variances for the suggested distribution on A have been calculated by using 
the programme DGQUAD from the library CERNLIB [12| and the results 



are presented in Table This Table shows that the series of properties of 
Poisson distribution (E'^ = A, Z)^ = A) take place only when the value of 
the parameter A > 3. 

It is appropriate at this point to say that 

oo POD A^6~'^ 

fix, X)dx = / dx = e-\{\). (4.3) 

Jo T[x + l) 

The function 

poo \x 

Jo r(x + i) 



is well known and, according to ref. ||l 

oo ■> n 

E ^7— Yy + 0(|Ar^-°-^) = e^ + 0(|Ar^) (4.5) 

n=—N ^ ' 

if A — > oo, \arg\\ < ^ for any integer N. Nevertheless we have to use 
the function 1{X) in our calculations in Sect. ^ and Sect. ^. We consider it 
as a mathematical trick to illustrate a possibility of constructing confidence 
intervals numerically for the real value x. 

Another approaches are also possible. At first, if we introduce a prior 

gA roo 

g{X) = , then we have equality / g{\)f{x, \)dx = 1 by natural way. 
y{X) Jo 

Also we can numerically transform the function /(x. A) in the interval x G 
(0, 1) so that 



roo roo roo 

/ f{x,X)dx = l, E^= xf{x,X)dx = X, DC= {x - E^f f{x, X)dx = X 
Jo Jo Jo 

(4.6) 

for any A. In these cases we can construct the confidence interval without 
introducing 1{X). 

Let us construct the central confidence interval for the continued ana- 
logue of Poisson distribution using the function /(A). 
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Table 3: The function 1{X), mean and variance versus A. 



A 


^(A) 


mean (-c/^j 


variance [Vt^) 


0.001 


-0.297 


-0.138 


0.024 


0.002 


-0.314 


-0.137 


0.029 


0.005 


-0.340 


-0.130 


0.040 


0.010 


-0.363 


-0.120 


0.052 


0.020 


-0.388 


-0.100 


0.071 


0.050 


-0.427 


-0.051 


0.113 


0.100 


-0.461 


0.018 


0.170 


0.200 


-0.498 


0.142 


0.272 


0.300 


-0.522 


0.256 


0.369 


0.400 


-0.539 


0.365 


0.464 


0.500 


-0.553 


0.472 


0.559 


0.600 


-0.564 


0.577 


0.653 


0.700 


-0.574 


0.681 


0.748 


0.800 


-0.582 


0.785 


0.844 


0.900 


-0.590 


0.887 


0.939 


1.00 


-0.597 


0.989 


1.035 


1.50 


-0.622 


1.495 


1.521 


2.00 


-0.639 


1.998 


2.012 


2.50 


-0.650 


2.499 


2.506 


3.00 


-0.656 


3.000 


3.003 


3.50 


-0.656 


3.500 


3.501 


4.0U 


-0.647 


A Ann 
4.0U0 


0.999 


4.50 


-0.628 


4.500 


4.498 


5.00 


-0.593 


5.000 


4.997 


5.50 


-0.539 


5.500 


5.497 


6.00 


-0.466 


6.000 


5.996 


6.50 


-0.373 


6.500 


6.495 


7.00 


-0.262 


7.000 


6.995 


7.50 


-0.135 


7.500 


7.494 


8.00 


0.000 


8.000 


7.993 


8.50 


0.000 


8.500 


8.496 


9.00 


0.000 


9.000 


8.997 


9.50 


0.000 


9.500 


9.498 


10.0 


0.000 


10.00 


9.999 
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5 The Central Confidence Intervals for the Con- 
tinued Analogue of Poisson Distribution. 



As we have noticed, for the discrete Poisson distribution the mean coincides 
with the estimation of parameter value, i.e. X = x. This is not true for 
a small value of A in the considered case (see Table ||). That is why in 
order to find the estimation of A(x) for a small value x it is necessary to 
introduce the correction in accordance with Table ^ Let us construct the 
central confidence intervals using a conventional method assuming that 

f{xX)dx = ^ (5.1) 

for the lower bound Ai and 

r f{x:\2)dx=^ (5.2) 

for the upper bound A2 of the confidence interval. 

Fig. m shows the introduced distributions (Sect. ^) with parameters 
defined by the bounds of confidence interval (Ai = 1.638, A2 = 8.493) 
for i; = A = 4 and the Gamma distribution with parameters a = 1, 
x + l = x + l = 5. 

The bounds of confidence interval with a 90% confidence level for the pa- 
rameter of continued analogue of Poisson distribution for different observed 
values X (first column) were calculated and are given in the second column 
of Table |. 

As a result (Table Q) the suggested approach allows one to construct 
confidence intervals for any real and integer values of the observed number 
of events for the values of parameter A > 3. Table ^ illustrates that the 
left bound of central confidence intervals is not equal to zero for small x. It 
shows that in this case a central confidence interval is not suitable. 

To anticipate a little, note that 90% of the area of Gamma distributions 
with the parameter x + 1 = x + 1 are contained inside the constructed 90% 
confidence intervals for the observed value x. However, for small values of 
X we have got values of the area close to 88%, i.e. less than 90%. 

The main goal of the proposed construction is to demonstrate a possibil- 
ity of using a continuous two-dimensional function (4.1) for the construction 
of confidence intervals in a frequentist meaning. 
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Table 4: 90% C.L. intervals for the Poisson signal mean A for total events 
observed x. 





bounds 


(Section 5) 


bounds 


(Section 6) 


X 


Ai 




Ai 


A2 


0.000 


0.121E-08 


2.052 


0.0 


2.303 


0.001 


0.205E-08 


2.054 


0.0 


2.304 


0.002 


0.292E-08 


2.056 


0.0 


2.306 


0.005 


0.666E-08 


2.061 


0.0 


2.311 


0.02 


0.218E-06 


2.098 


0.0 


2.337 


0.05 


0.765E-05 


2.166 


1.66E-05 


2.389 


0.10 


0.137E-03 


2.275 


2.23E-05 


2.474 


0.20 


0.186E-02 


2.490 


6.65E-05 


2.642 


0.30 


0.696E-02 


2.692 


1.49E-04 


2.806 


0.40 


0.161E-01 


2.891 


2.60E-03 


2.969 


0.50 


0.295E-01 


3.084 


5.44E-03 


3.129 


0.60 


0.466E-01 


3.269 


1.35E-02 


3.290 


0.70 


0.673E-01 


3.450 


2.63E-02 


3.452 


0.80 


0.911E-01 


3.629 


4.04E-02 


3.611 


0.90 


0.1179 


3.804 


6.12E-02 


3.773 


1.0 


0.1473 


3.977 


8.49E-02 


3.933 


1.5 


0.3257 


4.800 


0.2391 


4.718 


2.0 


0.5429 


5.582 


0.4410 


5.479 


2.5 


0.7896 


6.340 


0.6760 


6.220 


3.0 


1.056 


7.076 


0.9284 


6.937 


3.5 


1.340 


7.792 


1.219 


7.660 


4.0 


1.638 


8.493 


1.511 


8.358 


4.5 


1.946 


9.188 


1.820 


9.050 


5.0 


2.264 


9.869 


2.120 


9.714 


5.5 


2.590 


10.55 


2.453 


10.39 


6.0 


2.924 


11.21 


2.775 


11.05 


6.5 


3.264 


11.87 


3.126 


11.72 


7.0 


3.609 


12.53 


3.473 


12.38 


7.5 


3.961 


13.18 


3.808 


13.01 


8.0 


4.316 


13.82 


4.160 


13.65 


8.5 


4.677 


14.46 


4.532 


14.30 


9.0 


5.041 


15.10 


4.905 


14.95 


9.5 


5.406 


15.73 


5.252 


15.56 


10. 


5.779 


16.36 


5.640 


16.21 


20. 


13.(35 


28.49 


13.50 


28.33 
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6 Confidence Intervals for the Parameter of Pois- 
son Distribution in case of the real value of ob- 
served number of events. 

As follows from formulae (2.7) and (2.8) (see Figs. ^^) the probability of 
true value of parameter of Poisson distribution to be A in case of observed 
integer value x > distributes in accordance with the Gamma distribution 
with the parameters a = 1 and x + 1 = x + 1, i.e. according to formula (2.4) 

P{X\x) = ^e-\ 
x\ 

The possibility of constructing the continued analogue of Poisson distribu- 
tion suggests to assume that the Gamma distribution of true value of the 
parameter A takes place in case of the real value x > too (Figs. 
This supposition allows one to choose a confidence interval (for example) 
of a minimum length of all the possible confidence intervals of the given 
confidence level. The bounds of minimum length area, containing 90% of 
the corresponding area of probability density of Gamma distribution, were 
found numerically for several values of x. We took into account that x = x 
and required that Ai > 0. The results are presented in the third column of 
Table |. 

7 Conclusion 

In the paper the frequentist approach to construct the confidence interval for 
Poisson distribution parameter is considered. It is shown that the formula 
P(Ai < A < \2\x) = P{n < x\Xi) — P{n < x\\2) is a self-consistent defi- 
nition of the confidence interval in this case. It means that the probability 
of true value of parameter of Poisson distribution to be equal to the value 
of A in the case of m measurements rri, X2, . . . , Xm has probability density of 
Gamma distribution F^ 1+^" Xi ■ '^^^ results of constructing the frequen- 
tist confidence intervals for the parameter A^ of Poisson distribution for the 
signal in the presence of background with the known value of parameter A;, 
are presented. It is shown that the used procedure has both the Bayesian 
and frequentist interpretations. Also the attempt of introducing a continued 
analogue of Poisson distribution for the construction of confidence intervals 
for the parameter A of Poisson distribution is discussed. Two approaches 
are considered. Confidence intervals for different integer and real values of 
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the number of the observed events for the Poisson process in the experiment 
with a given integral luminosity are constructed. 
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Figure 1: The behaviour of the probabihty density of the true value of parameter 

A for the Poisson distribution in case of x observed events versus A and x. Here 
A^ 

/(x, A) = — -e^^ is both the Poisson distribution with the parameter A along the 
x' 

axis X and the Gamma distribution with a shape parameter a; + 1 and a scale 
parameter 1 along the axis A. 
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Figure 2: The Poisson distributions f{x, A) for A's determined by the confidence 
hmits Ai — 1.51 and A2 — 8.36 in case of the observed number of events x — A 
are shown. The probabihty density of Gamma distribution with a scale parameter 
a = 1 and a shape parameter a; + l = i; + l = 5is shown within this confidence 
interval. 
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Figure 4: Two-dimensional representation of the function f{x,X) versus A and x 
for the values f{x, A) < 1. 
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Figure 5: The probability densities /(x, A) of continued analogue Poisson distribu- 
tion for A's determined by the confidence limits Ai and A2 in case of the observed 
number of events x = 4 and the probability density of Gamma distribution with 
parameters a = 1 and x + l = x + l = 5. 
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